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1.1 Introduction 

The exciting physics of quantum phase transitions has been explored extensively in 
the last few years [1 Tl. The non-equilibrium dynamics of a quantum system when 
quenched very fast 13] or slowly across a quantum critical point [4 , 5 1 has attracted 
the attention of several groups recently. The possibility of experimental realizations 
of quantum dynamics in spin-1 Bose condensates [6 | and atoms trapped in optical 
lattices |i7j|8j| has led to an upsurge in studies of related theoretical models |l3]|4]|5]|9] 
[l0l[n][l2[l3][Il[l2[I3[ni[ISIl9l|20l|2T]|2^ 

In this review, we concentrate on the dynamics of quantum spin chains swept 
across a quantum critical or multicritical point or along a gapless line by a slow 
(adiabatic) variation of a parameter appearing in the Hamiltonian of the system. Our 
aim is to find the scaling form of the density of defects (which, in our case, is the 
density of wrongly oriented spins) in the final state which is reached after the system 
is prepared in an initial ground state and then slowly quenched through a quantum 
critical point. The dynamics in the vicinity of a quantum critical point is necessarily 
non-adiabatic due to the divergence of the relaxation time of the underlying quantum 
system which forces the system to be infinitely sluggish; thus the system fails to 
respond to a change in a parameter of the Hamiltonian no matter how slow that rate 
of change may be! 

We first recall the Kibble-Zurek argument ll34l l35l which predicts a scaling form 
for the defect density following a slow quench through a quantum critical point. We 
assume that a parameter g of the Hamiltonian is varied in a linear fashion such that 
g ~ gc ^ t /t, where g — gc denotes the value of g at the quantum critical point and 
T is the quenching time. Our interest is in the adiabatic limit, T ^ oo. The energy 
gap of the quantum Hamiltonian vanishes at the critical point as {g — gcY' whereas 
the relaxation time ^t, which is inverse of the gap, diverges at the critical point. It is 
clear that non-adiabaticity becomes important at a time i when the characteristic time 
scale of the quantum system (i.e., the relaxation time) is of the order of the inverse 
of the rate of change of the Hamiltonian; this yields 
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f^^.(f)->f^fe-g,r'-~ (^-j . (1.1) 

This leads to a characteristic length scale ^ given by 

(1.2) 

where £, is the correlation length or healing length. The healing length typically 
denotes the length over which a single defect is present. The density of defects in a 
(/-dimensional system scales as l/(^'' which leads to the Kibble -Zurek scaling form 
for the density of defects n given by 

„^T-^v/(vz+l)^ (1.3) 

The interesting aspect of the Kibble-Zurek prediction is that the scaling form of the 
defect density in the final state of a driven quantum system varies in a power-law 
fashion with the rate of quenching 1 /t, and the exponent of the power-law depends 
on the spatial dimension d and the static quantum critical exponents v and z- 

The Kibble-Zurek scaling has been verified in various exactly solvable spin mod- 
els and systems of interacting bosons ||4][T3][l6l[l7][TTl|51; it has been generalized 
to quenching through a multicritical point ||3T1 . across a gapless phase |20, 23J, and 
along a gapless line 1261 [28l . and to systems with quenched disorder ||T4| . white 
noisepTSl, infinite-range interactions fTl], and edge states |30|. Studies have also 
been made to estimate the defect density for quenching with a non-linear form 12111 . 
an oscillatory variation of an applied magnetic field 1321 or under a reversal of the 
magnetic field |33|. In the article in this book by Mondal, Sengupta and Sen (3E\, the 
quenching dynamics through a gapless phase and quenching with a power-law form 
of the change of a parameter as well as the possibility of experimental realizations 
have been discussed in detail. It should be mentioned that in addition to the density 
of defects in the final state, the degree of non-adiabaticity can also be quantified by 
looking at various quantities like residual energy 151 [37l and fidelity ID. 



1.2 A spin model: transverse and anisotropic quenching 

In this article, we will focus mainly on the one-dimensional anisotropic XY spin- 1/2 
chain in a transverse field l38l[39l which is represented by the Hamiltonian 



^ = - I [AcT>^+i + Jyaya:^, + ha;,i (1.4) 

n 

where a's are the usual Pauli matrices. (We will set Planck's constant h = 1). The 
spectrum of this Hamiltonian can be found exactly by first mapping the a matrices 
to Jordan-Wigner fermions c„ as 



n-l 

c„ = exp(/7r^ (7+(J^^), 

a~ = 2clc„~l, (1.5) 
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where — O'^ ± ial are spin raising and lowering operators respectively. These 
Jordan- Wigner fermion operators follow the usual anticommutation rules 

{c,;, , C„ } = 5,nn ; {c,„,C„} = Q ^ {cl„,cl}. (1.6) 

Applying a Fourier transformation to the Jordan- Wigner fermions along with the 
periodic boundary conditions, the Hamiltonian in ( 11.41 ) can be rewritten as 

H=-Y,{ [(7v + A)cosfe + /i] (4c^ + cl^c_^) 

*:>0 

+ /(/i-/v)sinfc {clc^-c^kCk}, (1.7) 

where k lies in the range [0, n\. The Hamiltonian is quadratic in fermion operators 
and hence exactly solvable using an appropriate Bogoliubov transformation ||39l . 
Using the basis vectors |0) (where no fermions are present) and \k^~k) — c^c^^jO), 
we can recast the Hamiltonian for wave number in a 2 x 2 matrix form 

[-i{J^-Jy)iimk ~h-{J^ + Jy)coiik\' ^ ' 

The spectrum for this is given by 

= [/!2+y?+7-,2 + 2/!(7,- + y,.)cosA: + 2/,y,.cos2A:]'/2. (1.9) 

This Hamiltonian has a very rich phase diagram as shown in Fig. II. II The vanishing 
of the energy gap for h = ±{Jx + Jy) for the critical modes at k = n and 0, respec- 
tively, signals quantum phase transitions from a ferromagnetically ordered phase to 
a paramagnetic phase with critical exponents v = z~ 1 [38 1. On the other hand, the 
vanishing energy gap at = Jy and h < Jx+ A denotes the anisotropic transition 
which marks the boundary between the two ferromagnetic phases with the critical 
exponents v and z being identical to the Ising transition. The meeting points of these 
two transition lines at h = ±{Jx + Jy) and Jx — Jy are multicritical points. 
Let us initiate our discussions with two types of quenching schemes: 

(i) quenching the magnetic field hast/T which we call transverse quenching ifTTlfTsjj . 
and 

(ii) the quenching of the interaction Jx as f /t which is referred to as anisotropic 
quenching 1 16|. 

In the process of anisotropic quenching, the system can be made to cross the multi- 
critical points A and B shown in the phase diagram. The quenching through a multi- 
critical point is a cardinal point of our discussion. We shall also discuss the quenching 
scheme where the anisotropy parameter y — Jx — Jy is quenched in a linear fashion 
keeping the system always on the gapless line h ~ Jx+Jy It will be shown that in 
either case, we arrive at new scaling behaviors for the defect density which cannot be 
obtained by a simple fine tuning of the Kibble-Zurek scaling form n ^ 1/ 

The reduction of the general Hamiltonian to a direct product of 2 x 2 matrices 
as given in (11.81 ) facilitates the application of the Landau-Zener (LZ) transition for- 
mula ll40l to estimate the non-adiabatic transition probability on passing through a 
quantum critical point. The general LZ Hamiltonian given by 
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Fig. 1.1. The phase diagram of the anisotropic XY model in a transverse field in the hjiJx + 
Jy) — y plane, where y = {Jx — J\) / {Jx+ Jy ) ■ The vertical bold lines given by /; / {Jx +Jy) = 
±1 denote the Ising transitions. The system is also gapless on the horizontal bold line y = 
for \h\ < Jx + Jy. FMx (FMy) is a long-range ordered phase with ferromagnetic ordering in 
the X (y) direction. The thick dashed line marks the boundary between the commensurate 
and incommensurate ferromagnetic phases. The thin dotted lines indicate the adiabatic and 
impulse regions when the field h is quenched from — oo to oo. The two points with coordinates 
y = and h/ (Jx + Jy) = ±1 denoted by A and B are the multicritical points. 



H = ei|l)(l| + £2|2)(2| + Zi(|l)(2| + |2)(1|) (1.10) 

closely resembles the reduced 2x2 spin Hamiltonian given in ( 11.8b . In the LZ Hamil- 
tonian, we set ei — £2 = f /t. The two energy levels zty/e^ + A^, where £1 = — £2 = £, 
approach each other with a minimum gap 24 at f = as shown in Fig. II. 21 The sys- 
tem is prepared in its initial ground state 1 1 ) at time f — 00 and should reach the final 
ground state |2) at f — > +00 if the dynamics is adiabatic throughout. A general state 
during the time evolution can be written in the form | i^(r)) = Ci (/) 1 1 ) + C2 (f ) |2) with 
the initial condition \Ci{t ~°°)|^ = 1- The LZ probability for the non-adiabatic 
transition is given by 1401 [37l 

= |Ci(r ^ oo)|2 =e-2'^^Vl|(^I-^2)l. (1.11) 

It is to be noted that the above formula is valid for a linear variation of the bare levels 
£1 and £2 and also when the off-diagonal term does not include any time dependence. 

Let us now discuss the transverse quenching of the XY spin chain discussed 
above. The phase diagram in (1.1) shows that when the transverse field is quenched 
from —00 to +00, the system is swept across the two Ising transition Unes at h = 
±{Jx + Jy)- Referring to the Hamiltonian in ( 11.81 ). it is clear that when h is quenched 
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Fig. 1.2. The two levels correspond to the energy levels varying with time; 24 is the minimum 
gap between these two levels, and |1) and |2) are the eigenstates in the asymptotic limits 
t ±oo. The ground state changes its characteristics from 1 1) at / ^ — oo to |2) at f ^ +<». 



as /i ^ f/r ifTSll . the dynamics of the spin chain effectively reduces to a LZ problem 
in the two-dimensional reduced Hilbert space spanned by the basis vectors |0) and 
\k,—k), with the initial ground state being |0). Here, the diagonal and off-diagonal 
terms are denoted by = ±[h + {Jx + 7,. ) cos A:] and A — i {Jx — /y ) sin A; respectively. 
Denoting a general state vector \lfk{t) at an instant t as \'^k{t)) = Ci jt|0) +C2.a.|A:, —k) 
with \C\^k[t ~> " 1' the non-adiabatic transition probability for the mode k is 

directly obtained using the LZ transition formula 

Pk = \Cx.k{t - = (1.12) 

The density of defects (i.e., density of wrongly oriented spins) in the final state is 
obtained by integrating pi^ over the Brillouin zone ifTSll . 




(1.13) 



In the adiabatic limit x ^ the transition probability is non-zero only for modes 
close to the critical modes k = Q and tt; this allows us to extend the limits of inte- 
gration from — oo to +00. Noting that the critical exponents v = z = 1 for the Ising 
transition and d —\ \ hence the 1 / -^/t scaling of the defect density is consistent with 
the Kibble-Zurek prediction. 

On the other hand, if we look at the anisotropic quenching Jx ^ t/T, the off- 
diagonal terms of the Hamiltonian in dl.SI ) pick up a time dependence and hence a 
direct application of the LZ transition probability is not possible. To overcome this 
problem, we rewrite the Hamiltonian in the basis of the initial and final eigenstates 
when Jx —°° and Jx °°. The eigenstates of the Hamiltonian in these limits are 
given by 

|eu> = sin(V2)|0) +/cos(V2)|^,-^) 
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and 



\e2k) =cos(A:/2)|0) -/sin(A:/2)|A:,-fe), 



with eigenvalues Ai —t/x and A2 = — f/t respectively; the system is in the state \e\k) 
initially. A general state vector can be expressed as a linear combination of jei^.) and 
\e2k), 

\Wk{t))^Cydt)\eik)^C2k{t)\e2k)- (1.14) 

The initial condition in the anisotropic case is C\k{—°°) = 1 and C2a(— °°) = 0. The 
unitary transformation to rewrite the Hamiltonian in the \en) and \e2k) basis is given 
by Hj^{t) = U'^Hk{t)U, where 



U 



and the new Hamiltonian is 



cos(A;/2) sin(A:/2) 
— isin(A;/2) /cos(A;/2) 



H[{t) = -[/!+ (y^+7y)cosfe]/2 

Jx + JyCos2k + hcosk /v sin 2A: + /1 sin A: 
Jy s,m2k + hs,mk —Jx — Jy cos2fc — hcosk 



(1.15) 



By virtue of the unitary transformation, the time dependence is entirely shifted to 
the diagonal terms which makes it possible to apply the LZ transition formula. Eval- 
uating the probability of a non-adiabatic transition for the mode k and integrating 
over the Brillouin zone, it can be shown that the density of defects in the final state 
following the anisotropic quenching is given by |fT6l 

47, 



h 

for h>2Jy. 



(1.16) 



The scaling behavior of the density of defects with T is shown in Fig. II. 31 

Although Eqs. (11.16b satisfy the Kibble-Zurek scaling, the density of defects di- 
verges for h ~ 2Jy, i.e., in a passage through the multicritical points. This necessitates 
a generalization of the Kibble-Zurek prediction for quenching through a multicritical 
point. In the next section, we propose a generalized scaling for the density of defects 
valid for quenching through a critical point as well as a multicritical point. 

In Fig. 11.41 we present the variation of the von Neumann entropy density of the 
final state defined by |fT3l 



■s = - / — [PkMPk) + (1 - W)ln(l - w) ]. 

Jo 7Z 



(1.17) 



following the anisotropic quenching with the rate of quenching. Even though the final 
state is a pure state as a result of a unitary dynamics, it can also be viewed locally, or 
on a coarse-grained wave vector scale, as a decohered (mixed) state 1 1 3J . In the limit 
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Fig. 1.3. Plot of kink density n versus T as obtained for h = 0.2,0.4,0.6,0.8 (from bottom to 
top in the large T region), with Jy = \ for anisotropic quenching. For large T, n increases with 
increasing h, whereas for small T, it decreases with increasing h. 



of T ^ 0, the system does not get enough time to evolve; hence it largely retains its 
initial antiferromagnetic order which results in a low local entropy density. On the 
other hand, for slow quenching (t — > oo), the system evolves adiabatically, always 
remaining close to its instantaneous ground state; this results in a final state with 
local ferromagnetic ordering, and hence again with a low local entropy density. We 
observe that the entropy density shows a maximum at a characteristic time scale Tq 
where the magnetic ordering of the final state also changes from antiferromagnetic 
to ferromagnetic. 




Fig. 1.4. Plot of von Neumann entropy density s, staggered magnetization msx and magnetiza- 
tion m-x as a function of t/tq, for Jy = \ and h = 0.2 in the anisotropic quenching case. 
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Let us briefly mention a few other interesting variants of the transverse quench- 
ing scheme. In a repeated quenching dynamics of the same model, the transverse 
field h{=t It) is quenched repeatedly between —oo and +00 f25 1 . We refer to a sin- 
gle passage from h —oa to h — > +00 or the other way around as a half-period of 
quenching. An even number of half-periods corresponds to the transverse field being 
brought back to its initial value of —00, whereas, in the case of an odd number of 
half-periods, the dynamics is stopped at /z ^ +00. The probability of a non-adiabatic 
transition at the end of I half-periods can be shown to follow the recursion relation 

p,{l) - (1 - e-^-y) - (1 -2e-2-r) [1 _^,(; _ 1))] , (1.18) 
eventually yielding the simplified form of 



For large t, the density of defects is generally found to vary as 1 / -/r for any 
number of half-periods. On the other hand, for small T, it shows an increase in 
kink density for even values of /. However, the magnitude is found to depend on 
the number of half-periods of quenching. For two successive half-periods, the defect 
density is found to decrease in comparison to a single half-period, suggesting the ex- 
istence of a corrective mechanism in the reverse path. The entropy density increases 
monotonously with the number of half-periods, and shows qualitatively the same be- 
havior for any number of half-periods. For a large number of repetitions (/ ^oa), the 
defect density saturates to the value 1/2 while the local entropy density saturates to 
In 2. 




0.001 I — • ^ — • — ^ — • — ^ — • ^ — • — ^ — • ^ 

Fig. 1.5. Plot of kink density n after / half cycles as a function of T, for Jx — Jy = i and 
/ = 1,2, 3,4. «(4) denotes the defect density as obtained from the analytical expression. In the 
limit of large -r, the numerical results match perfectly with the analytical results. 



The effects of interference on the quenching dynamics of Hamiltonian in (11.8b 
when the transverse field h{t) varies sinusoidally with time as h = hocoscot, with 
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|f I < n/(£) has also been studied in a recent work ll32l . In this scheme of quench- 
ing, the time interval between two successive passages through the quantum critical 
points can be small enough for the presence of non-trivial effects of interference in 
the dynamics of the system. It has been shown that for a single passage through a 
quantum critical point, the interference effects do not contribute leading to a kink 
density which goes as n ~ On the other hand, repeated passages through the 
quantum critical points result in an oscillatory behavior of the kink density as well 
as the entropy density. 

1.3 Quenching through a multicritical point 

We shall now propose a general scaling scheme valid for quenching through a mul- 
ticritical point as well as a critical point f3T| using the LZ non-adiabatic transition 
probability I40. ,37J discussed before. We begin with a generic (/-dimensional model 
Hamiltonian of the form 



H = V^'^(k) [(A(f) + fe(k)) a~ + zi(k) a+ + 4*(k) o'] v/(k), (1.20) 



where a* = ± io-' , b{^) and A (k) are model dependent functions, and i^(k) 
denotes the fermionic operators (y/i (k), v/2(k)). The above Hamiltonian represents, 
for example, a one-dimensional transverse Ising or XY spin chain f39l as shown in 
Eq. (1.8), or an extended Kitaev model in d — 2 written in terms of Jordan- Wigner 
fermions 14111201 . The excitation spectrum takes the form 



We assume that the parameter X{t) varies linearly as f/r and vanishes at f = 0. 
The parameters b{Vi.) and A (k) are assumed to vanish at the quantum critical point in 
a power-law fashion given by 



where we have taken the critical mode to be at ko ~ without any loss of generality. 
Eq. ( II.22I 1 also implies that the system crosses the gapless point at f when X 
and biVi) and 4(k) also vanish for the critical mode k = 0. Many of the models 
described by Eq. ( 11.20b exhibit a quantum phase transition with the exponents asso- 
ciated with the quantum critical point being z\> zi and hence z — zi — 1. We shall 
however explore the more general case encountered at a multicritical point where the 
dynamically exponent z is not necessarily given by 22- 

The Schrodinger equation describing the time evolution of the system when X 
is quenched is given by /(9\//k/(9f ~ Hy/^, where we have once again defined i/k as 

= Ci,k|0) +C'2.k|k,— k). Using the Hamiltonian in Eq. (11.201 1. we can write 



k 




(1.21) 



b{k) - |k|'i and A{k) - |k|^2^ 



(1.22) 




.<5C2,k 



( 




(1.23) 



dt 
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One can now remove b{k) from the above equations by redefining f/T + /7(k) —tt/x; 
thus the exponent z\ defined in Eq. ( 11.22b does not play any role in the following 
calculations. Defining a new set of variables Ci k = Ci,kexp(// dt' t' jx) and C2,k = 
C2,kexp(— ! /' dt' t' I t), we arrive at a time evolution equation for C\ (k) given by 

Further rescaling f tx^l^ leads to 

If the system is prepared in its ground state at the beginning of the quenching, i.e., 
C\ (k) = 1 at f = — oo, the above equation suggests that the probability of the non- 
adiabatic transition, pk ~ lim(^+„o |Ci kp, must have a functional dependence on 
|4(k)pT of the form 

= /(|4(k)|2T). (1.26) 

The analytical form of the function / is given by the general LZ formula ll37l l40l . 
The defect density in the final state is therefore given by I 31J 




(1.27) 



The scaling |k| jkp^^-j- finally leads to a scaling of the defect density given by 

n - 1/t''/(2^'2). (1.28) 

Let us recall the quenching dynamics of the transverse XY spin chain discussed 
before lfT3l [TFl . When the system is quenched across the Ising or anisotropic crit- 
ical line by Unearly changing /z or 7^ as f /t, 4 (k) vanishes at the critical point as 
A (k) '--^ |k| yielding Z2=z=\\ hence the generalized scaling form given in Eq. ( 11.28b 
matches with the Kibble-Zurek prediction with v z = 1 as mentioned before. The 
situation, however, is different when the system is swept across a multicritical point. 

Putting h — 2Jy in the Hamiltonian in ( 11.151 ), let us analyze the scaling of the 
diagonal and off-diagonal terms of the Hamiltonian at the multicritical point. When 
Jx — Jy, the diagonal term of the Hamiltonian scales as —Jy{n — k)'^ whereas the off- 
diagonal term \A{k)\ ~ \n — k\^ when expanded around the critical mode at k — n. 
The dynamical exponent is obtained from the diagonal term so that z = Zi — 2. As 
discussed above, the off-diagonal term or more precisely the exponent Z2 determines 
the scaling of the defect density. The exponent Z2 = 3 for the XY multicritical point; 
hence the defect density scales as 1 / T^/^. The density of defects obtained by numer- 
ical integration of the Schrodinger equation in ( 11.23b as shown in Fig. 1 1.61 supports 
this scaling behavior. For a non-linear quench across a multicritical point [36], when 
the parameter A is quenched as A ^ {t /x)"'sgn{t), both z\ and Z2 come into play 
and the scahng of defect density gets altered to n ^-dav I[a(z2v+i)+ziv(\-a)].^ ^yAs 
reduces to the form presented above for a~ \ and ziV = 1. 
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Fig. 1.6. n vs T obtained by numerically solving Eq. ( I1.23t at the multicritical point of the XY 
model in a transverse field, with h = 10 and Jy = 5. The line has a slope of —0. 16. 



1.4 Quenching along a gapless line 

We now shift our focus to quenching the XY Hamiltonian in il.7\ along the gapless 
line h = +Jy by varying the anisotropy parameter Y—Jx — A as f /t, keeping +Jy 
fixed ||26| . For convenience, let us set Jx + Jy = h — 1. The excitation gap vanishes 
along this line for the mode k = n. We rewrite Eq. (ll.SI l with the present notation in 
the form 

1+cosA: iysAwk 
—iysink — 1— cosA: 



H, 



(1.29) 



We once again encounter a situation in which the off-diagonal terms are time- 
dependent. Noting that the asymptotic form of the Hamiltonian at f ^ ±oa is given 

by 

H=-sink&-\ (1.30) 

T 

we make a basis transformation to a representation in which a" is diagonal. The 
Hamiltonian in ( |1.29t then gets modified to the form 



(f/r) sinA: 
1 +cosk 



1 +cosA: 

- (f/r) sin A: 



(1.31) 



where the time dependence has been shifted to the diagonal terms only. Applying the 
LZ transition probability formula, the probability of excitations is found to be 



Pk ■■ 



-nx\Ai^\~ / sink 



(1.32) 



where = 1 1 +cosA:p — {% — when expanded about k= n. Integrating pi^ 
over the Brillouin zone, we find the that density of defects falls off as « ^ 
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Fig. [LT] which shows n vs t obtained by numerical integration of the Schrodinger 
equation confirms the t^^/^ behavior. 

It is to be noted that the Kibble-Zurek formalism cannot address defect genera- 
tion along a gapless line. We will now present another general scaUng argument for 
moving along a gapless hne in a (/-dimensional system. 




log X 

Fig. 1.7. nvs T obtained by numerically solving Schrodinger equation with Y(t) ^t/t keeping 
h = 2Jy. Also shown is a line with slope 1/3 for comparison. 



Let the excitations on the gapless quantum critical line be of the form Ok ^ A |k|^, 
where z is the dynamical exponent and the parameter A —t/iis quenched from — oo 
to °o. Using a perturbative method involving the Fermi Golden rule along with the 
fact that the system is initially prepared in the ground state, the defect density can be 
approximated as lil9J 



OA 



(1.33) 



Let us assume a general scaling form for the instantaneous excitation, 5(Bk(A') = 
k"f{^), where k — |k| and k" denotes the higher order term in the excitation spec- 



trum on the gapless line. Defining a new variable ^ — Xk^ 
behavior of the defect density as 1 26 1 



we obtain the scaling 



(1.34) 



The case d — I, a — 2 and z = 1 has been discussed in this section. Note that the 
correlation length exponent v does not appear in the expression in Eq. (11.34b because 
our quench dynamics always keeps the system on a critical line. The scaling for 
quenching along a gapless line is less universal in comparison to the general Kibble- 
Zurek prediction. For a non-linear quench |22| A ~ {t /T)"sgn{t), this scaling form 
gets modified to n ~ '^-da/la{a+i)-z\.^ jj^jg reduces to l/x^^^ foi a = d = z = I and 
a = 2. 
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In passing, let us introduce a variant of the XY Hamiltonian with an alternating 
field given bv Il42ll28l 



H = -]-[Y^(J, + Jy){a]a]^, 
^ j 

+ (J:,~Jy){a]a]+,~a]<y]^^ + ih-{-iy5)aj]. (1.35) 

The strength of the transverse field alternates between h + 5 and /i — 5 on the odd 
and even sites respectively. For 5 = 0, we recover the conventional XY model in a 
transverse field. This Hamiltonian can also be solved exactly by a Jordan-Wigner 
transformation by taking care of the two underlying sublattices, i.e., by defining 
Uk and bk as two different Jordan-Wigner fermions. The Hamiltonian in the basis 
_k) can be written as 



{a'k,a-k,bl,b- 



H, 



h + Jcosk iysink —5 

—iysink —h — Jcosk 5 
5 Jcosk — h iysink 

— d —iysink —Jcosk + li 



(1.36) 



The excitation spectrum of the above Hamiltonian is 

A± ^ [h^ + d^+J^cos^k + fsiir-k 



±2 J/1252 + /1272 cos2 ;t + 5272 sin2 k ] 



1/2 



(1.37) 



where J = Jx + Jy and 7 = 7, — Jy, with four eigenvalues given by . In the ground 
state, —Aj^ and — A^r are filled. The vanishing of A^r dictates the quantum critical 
point, and the critical exponents are obtained by studying the behavior of Aj^ in the 
vicinity of the critical point. The minimum energy gap in the excitation spectrum 

52+72 and 



occurs at A; = and k = n/2. The corresponding phase boundaries 

= + 72 signal quantum phase transitions from a paramagnetic to a ferromag- 
netic phase respectively. We extend the study of quenching through a gapless phase 
in this Hamiltonian by varying 7 as before along the phase boundary = 8^ +J^. 
On this gapless line, the dispersion of the low-energy excitations at ^ can be 
approximated as 



A,: 



j^k^ 



y2j2l^2 



(1.38) 



Y 4(52+72) 52+y2 ■ 
This suggests a truncation of the Hamiltonian in ( |1.36l l to a 2 x 2 LZ Hamiltonian 



f{t)k Pk^/2 
Pk^/2 -Y{t)k 



(1.39) 



where 7 and 7 are renormalized parameters given by 7 = 77/ ^/5^ and 72 ~ 
72/ V 52 +72. The argument which justifies the truncation of a 4 x 4 matrix to a 
2x2 matrix is presented in Ref. [26 1. The diagonal terms in Eq. (|1.39l l describe 
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two time-dependent levels approaching each other linearly in time (since 7 = f /t), 
while the minimum gap is given by the off-diagonal term J^k^ /2. The probability of 
excitations pi^ from the ground state to the excited state for the A;-th mode is given by 
the LZ transition formula 1,40, 37] 

and we get back the t^'/^ scaling form of the defect density as discussed above. 



1.5 Conclusions 

The non-equilibrium dynamics of a quantum system driven through a quantum crit- 
ical point is indeed an exciting and fascinating area of research. The possibility of 
experimental realizations of quenching dynamics adds to the importance of the theo- 
retical research. In this review we have discussed the scaling relations for the defect 
density in the final state following a slow quench across a multicritical point and 
along a gapless line. In both cases, the relations obtained here using the LZ transi- 
tion formula are not directly derivable using the Kibble-Zurek argument. But in all 
the cases, the defect density scales as a power law with the rate of quenching. It 
should be noted that the idea of a dominant critical point has also been invoked to 
justify the scaling along a gapless line in one-dimensional spin models [28|. A gen- 
eralization of the studies presented here to higher dimensional quantum systems and 
to systems with quenched disorder are some of the future directions of this subject. 
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